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We propose a new holonomy formulation for spin foams, which naturally 
extends the theory space of lattice gauge theories. This allows current 
spin foam models to be defined on arbitrary two-complexes as well as 
O ■ to generalize current spin foam models to arbitrary, in particular finite 

groups. The similarity with standard lattice gauge theories allows to apply 
Ch^. standard coarse graining methods , which for finite groups can now be 

easily considered numerically. We will summarize other holonomy and spin 
network formulations of spin foams and group field theories and explain 
how the different representations arise through variable transformations in 
the partition function. A companion paper will provide a description of 
boundary Hilbcrt spaces as well as a canonical dynamic encoded in transfer 
operators. 

1. Introduction 

Spin foam models arc non-perturbativc candidate theories for quantum gravity [TJ [2] , 
which succeed in combining discreteness of spectra of geometric operators [3] [4] with 
the correct semi-classical behaviour of its elementary building blocks [5] [51 H] ■ One of 
the key challenges is to obtain a better grasp on the continuum limit of spin foams, 
obtained by considering the limit of (infinitely) many building blocks. This will not 
only answer questions on the large scale behaviour (which here we will understand as 
the one involving many building blocks) but also address the status of diffeomorphism 
symmetry in these models [HI [HI QZ3 [HI HI] • 

To investigate this kind of limit involving many building blocks rcnormalization and 
coarse graining methods seem to be most appropriate [131 HH US] US] . To this end 
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an understanding of the space of models, on which the renormalization flow takes 
place, is crucial. This allows also also to distinguish between relevant and irrelevant 
choices in the construction of spin foam models and helps to reveal the basic dynamical 
mechanisms of these models. 

In this work we therefore propose a unified description of spin foam models based 
on group or holonomy variables. In this form spin foam models appear to be similar 
to lattice gauge theories, have however a richer structure of dynamical ingredients. In 
the description with the fewest parameters our models arc defined by an edge function 
E and a face weight u>, whereas lattice gauge theories just feature the face weights 
lu. All the current spin foam models, namely Barrett-Crane, EPRL and FK models 
[HI UHl [H] , can be described by this choice. 

For this description to be useful for coarse graining one needs to show that these 
structures are stable under coarse graining (within some approximation scheme). In- 
deed we will investigate different coarse graining schemes for two and three dimensional 
models analytically and numerically and show that this space of models presented here 
supports a rich structure for the renormalization flow. To be able to apply numerical 
techniques we consider spin foam models with finite groups, a concept introduced in 
[20l ITS] . A further advantage of the holonomy description here is that models, anal- 
ogous to the Barrett Crane or EPRL model, can be easily introduced also for finite 
groups. This will ease the investigation of (for instance symmetry) properties of these 
models as well as the numerical investigation of the large scale limit. 

There arc different holonomy representations proposed already in the literature |211 
l22l l23l 124] , the relation to these will be explained in section [4] The advantage of 
the formulation presented in this work is that it features - compared to the other 
formulations - a minimal set of model parameters and makes maximal use of the 
common properties of the spin foam models proposed so far in the literature. Most 
prominently that the imposition of the simplicity constraints factorizes in a specific 
sense, see also [2]. 

We will also explain the relation to the so-called operator spin foam models [25], 
which share with the models proposed here, that they are defined on arbitrary two- 
complexes. Indeed the structure of the holonomy representation worked out here makes 
the definition on a two-complex completely natural. 

Other advantages of the holonomy representations introduced here are the natural 
definition of boundary Hilbert spaces and a straightforward relation to projected spin 
networks [55] HZ] as well as allowing the investigation of the semi-classical limit for 
arbitrary complexes by exploring the distributional character of the integration kernels 
of the partition functions. These features will be explained in companion papers 

EE1C!]. 

In section [3] we give the first main result of this paper, a definition of spin foam 
models as generalized lattice gauge theories written entirely in terms of holonomics. 
Section 2] discusses the relation to operator spin foams and other holonomy formula- 
tions for spin foams, including group field theories. The next section gives a description 
of the Barrett-Crane and EPRL spin foam models in the holonomy formulation. These 
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insights are used to generalize the Barrett Crane and EPRL models to finite and other 
Lie groups. 

In section [5] we discuss coarse graining of these models, in particular in two and 
three dimensions. We close with an outlook and summary in sections 16.41 and [71 



Here we shortly motivate some ideas behind the holonomy representation for spin 
foams, which will be explained in more detail later on. 

Let us start with reviewing standard lattice gauge theory. Such a theory requires 
the choice of a gauge group G and face weights w, which are class functions on the 
group. These face weights encode all the dynamics of the given model, as the partition 
function is given by 



On the other hand the face weights for i.e. lattice Yang Mills theory arc defined using 
the background metric of the lattice. This should however be avoided if one wants 
to quantize gravity, where the metric itself is a variable. Indeed, the starting point 
for spin foam model construction is BF theory where the face weights are given by 
u! = Sg the delta function on the group. This makes the partition function in general 
divergent (even for compact groups but not for finite ones) . This divergence indicates 
another gauge invariance, known as translation symmetry. BF theory is a topological 
field theory, which coincides with the zero coupling or zero temperature limit of lattice 
Yang Mills theory. 

To construct a non-topological theory and nevertheless keep the trivial face weights, 
we have to introduce more structure. First let us just double the variables from one 
group element g e per edge to one group element g ve per edge-vertex pair (or half edge). 
We define g ev — (gve) -1 - Then the partition function (TTJ) reads 



and we can regain the form ([T]) by variable redefinitions g e = g ve g ev > and performing 
\E\ trivial integrations. Here and in the following we assume a normalized and left 
and right invariant group integration measure. 

To actually change the model we introduce another set of variables, namely a group 
element g e f per face-edge pair. These are inserted into the face holonomies so that 
uif is now evaluated on gj = g ve g e jg ev ' <jVe' In this way the alternative face 
holonomy Gf = g ve g e v' <?u'e' '•• can nevertheless take non-trivial values even if we 
choose ujf = Sq. Allowing all possible g e f would however again trivialize the model. 
(This time it corresponds to the high temperature fixed point of lattice Yang Mills, 



2. A short motivation of the models 




(i) 
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where utf = 1 is the constant function.) A non-trivial model can be obtained by 
weighing the inserted group elements g e f via a function E. The final model will be 
given by 




(3) 

As will be explained in the following the partition functions of all current spin foam 
models can be expressed in this way. The choice of the function E determines the 
model, the face weights are usually fixed to be Sq- 

The models can be rewritten in several ways, we will present different formulations 
in section 2J Here we will just note that one could integrate out the g e f variables in 
the formulation to obtain effective face weights 

w f(9ev9ve' j 9e'v'9v'e" 7 ' ' ' ) = 




Such a formulation (for the Barrett Crane model) appeared for instance in [211 [22] 
and is nearest to standard lattice gauge theory. See also [22 1201 EI] for a geometric 
interpretation of the holonomy variables and simplicity constraints involved there. 

Note that the face weights are not just class functions of a face holonomy, as is the 
case in standard lattice gauge theory, but will be functions of a certain number of 
group elements. 

3. Holonomy Formulation of Spin Foam Models 

Let us describe the models in more detail. These will be defined on an oriented 2- 
complcx C with vertices v, edges e and faces /. The gauge group will be denoted G, 
which we will assume to be finite or compact and equipped with a normalized left and 
right invariant group (Haar) integration measure. 

The variables in which the partition function will be expressed are group elements 
g e f assigned to pairs e C / and g ve for v C e with the understanding that g ev = g~^. 
One can identify them as holonomies of a G-connection along certain paths within 
faces / (see Figure H]) 

From these data we can define the holonomy G/ 

Gf = g V e9ef9ev'9v'e' "• 9 e Wv (5) 

With the help of Figure [TJ one can readily see that G / describes the holonomy as- 
sociated to a geometrically trivial path, i.e. the choice uj = Sg as a weight to this 
holonomy appears most natural and is also used in the EPRL-FK or BC model. For 
rcgularization purposes, this 8q is often replaced by a heat kernel, see e.g. [32] . We 
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V V > V v > 

Figure 1: Holonomy variables associated to half edges (g ve ) and wedges (g e f) of a face. 



will therefore leave the face weight w as a parameter of the theory. As in standard 
lattice gauge theory uj is required to be a class function of the group, furthermore we 
ask that oj(g) — is satisfied for all g £ G. 

The class of models we arc considering arc specified further by a subgroup G' C G. 
We will choose the edge function E to be invariant under conjugation with elements 
in G' . From (j4j one can see that E functions which are invariant under conjugation 
of the full group, i.e. class functions, will lead to effective face weights which can be 
considered to just depend on the (reduced) face holonomy g ev g ve ' g e 'v' • • • ■ This case 
will therefore lead to a model equivalent to a lattice gauge theory (with face weights u / 
that might explicitly depend on the face, more precisely on the number of edges in the 
face). In this sense the choice of G' defines an invariance property of the effective face 
weights ([!]). As for the face weights we will furthermore require from E the property 
E{g) = E{g- 1 ) for g G G. 

Finally the state sum is defined as 

Z= fdg ev dg ef l[uj(G f )l[E(g ef ) . (6) 

J f ef 

For the choice u> = 5q the partition function may be divergent, even for compact 
groups. The definition is extendable to non-compact groups, but in this case there 
is an ambiguity in the normalization of the Haar measure and additional divergences 
due to gauge orbits. However, for finite groups the delta distribution is just a regular 
functions, thus the above is well defined. 

3.1. Parametrization of models 

The space of such holonomy spin foam models is described by a choice of the edge 
functions E and the face weights ujf. In the following we will give the spin representa- 
tion of these entities as such a representation also offers a convenient parameterization 
of the space of models. 

The face weight ojf is a class function and as usual can be expanded as 
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= ^2 d ^(p) u}P x P (g), (7) 

p 

where Xp — tr p (D p (g)) is the character of the unitary irreducible representation (irrep) 
p with D p (g) the representation matrix of g. 

The edge functions can be expanded in terms of the irreducible unitary represen- 
tations of G and G' , which we denote p and k respectively. To this end we utilize 
injection maps I(p, k)d' 

I(p,k) d e Hom G ,{p 7 k), (8) 

where the d index labels a basis of Home 1 (p, k) the space of G' covariant function 
from the representation space of p, the dual of p, to the G' irrep k. 

The map I(p, k)d identifies the representation k, which in case there are more than 
one representation is furthermore labelled by d, in the reduction of the irrep p over the 
group G' . The index d is hence only non-trivial in cases where there are more than 
one k representations included in p. 

We assume that the basis is normalized (up to a phase, not relevant in the following) 

I(p,k)ll(p,k) d , = 8 dd ,l k . (9) 

This allows to expand E as 

E(h) = dim(p)< dd , tr p (D p (h) I(p, k) d l\p, k) d ,) . (10) 

p 7 k 

Hence for given groups G, G' the choice of E function can be encoded into the param- 
eters e p k dd , . The index d is a degeneracy index for the case where the G' reducible 
representation p contains more than one copy of the irrep k. For G = SU(2) x SU(2) 
and G" = SU{2)di a g this case does not occur, and we have only to deal with the set 

4- 

E encodes the analogue of the simplicity constraints, arising in the Plcbanski for- 
mulation of gravity, for the spin foam model at hand, sec for instance [2]. As in the 
classical Plcbanski formulation, where the simplicity constraints are enforced with La- 
grange multipliers, we had to introduce here additional variables, the wedge holonomies 

9ef- 

If the face weight is given by the delta function, the product of group elements 
around the face, G/, is flat. This is however not the usual holonomy around the face 
Gf = g ve g ev > : g v > - e ' • • • , but G/ interwoven with the g e f. If we force g e f = 1 by choosing 
E = Sq we have G/ = G/, and obtain a theory of flat connections, that is BF theory. 
Choosing a different E function relaxes the constraints on flatness, just as we expect 
from the simplicity constraints. 

3.2. Boundaries and Hilbert spaces 

We can also include boundaries into our models, which are represented as graphs T in 
the two-complex C. The set of edges of Y will be denoted by Y e C C e and the set of 
vertices r„ C C v . We write Y ev for the set of pairs v G e in r„ x Y e . 
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By dropping integration over the group elements associated to the pairs of r e „ we 
arrive at the partition function 



z r (o[g ev }= j (n<w) I n d -9- 1 (n^/)J (iw^J • ^ 

which can now be interpreted as an element of the Hilbert space L 2 (G' re "'). (Integra- 
bility in this general form is only guaranteed for finite groups). 

As the integrand of the partition function is invariant under g ev — ¥ h' e g e v9v , h e f — > 
h' e 1 h e fh' e for h' e eG'cG and g v £ G, the partition function can be considered to 
live in a smaller subspace of L 2 (G' rc "'), that is, 



ZeH r =L 2 (G^ /fGl r "l x G' 1 ^ 1 )) . (12) 



This Hilbert space also describes the effective face weights w' defined in (UJ. 

We call this space the universal boundary space for the holonomy spin foams, which 
is fixed by a choice of graph, and the two groups G,G'. This Hilbert space and its 
relation to projected spin networks |26[ I27j will be further explored in the companion 
paper [25] . 



4. Relation to other formalisms 

The form of the state sum ((5]) was chosen such that it mimics the main features 
of the EPRL (and other known) spin foam models, but has a minimal set of free 
parameters. This appears to be a very convenient choice, but it is by far not the only 
one. There have been other forms used in the literature, see for instance [2T1 1221 123"l 124) . 
which include weights for different geometric objects, which might appear under coarse 
graining. The aim of this section is to give a precise relation between the .E-function 
formalism used in this paper, and the other formulations. 



4.1. Operator Spin Foam models 

The operator spin foam models (OSFM), introduced in |25j . were originally formulated 
with the aim of allowing for arbitrary operators P e associated to edges e of C. The 
state sum Z is evaluated by contracting the P e at the vertices v of C according to the 
combinatorics of the 2-complcx. 

Starting with a 2-complex with irreps pf of G associated to faces /, an OSFM is 
specified by operators mapping the intertwiner space onto itself 

P e : Inv G (p /l (g> . . .pf n ) — ► InvcGo/! ® • • -P/J, (13) 

as well as face weights tu Pf . There is a dual formulation of this model [23], using 
holonomics h e f G G for all pairs (e C /). This holonomy can be thought of as starting 
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Figure 2: Holonomy variables in the dual description of operator spin foam models. 



at one vertex v, going to the middle of the face /, and going to another vertex v', so 
that e goes from v to v' in the orientation induced by /. The Fourier transform of the 
P e operators is given by functions C e 



Pk \ k J 



(14) 



where we assume for simplicity that the edge orientation agrees with all orientations 
from all faces meeting at it. Then, the state sum in terms of holonomies is given by 



^OSFM 



Jdhej ll[u;(G f ))l[C e (h efl ,...,h efn ) (15) 



The Gf is the ordered product of holonomies h e f around a face /. This forms the 
same path as ([5]). 

The i?-function formalism is a special case of the OSFM, in which the edge op- 
erator P e factorizes over the edges, and the invariance property O e (hi, . . . , h n ) = 
C e {gLh\gR, . . . ,gLh n gn), which follows from (|13[) and (|14[) . is explicitly enforced by 
projectors onto the intertwiner space. One can write the i?-function state sum in 
the form (|15[) by performing the variable transformation 

h e f = gvegef9ev'- (16) 

This removes the g ev - dependence of the face functions lo : and one can carry out the 
integration over the g ev , defining 

C e (h e f h e f n ) ■■= / dg ev dg ev > ] [ E(g ev h e f k g v > e ) (17) 

G2 fkDe 

for each edge. This results in (fT5|) . as one can readily check. 

Hence one can see that the i?-function state-sum © is a special case of the OSFM 
state-sum (|15p , in which the operator P e , which is a map from the invariant subspacc 
of pj x ® • • ■ pf n to itself, can be written as a part which factorizes over every face, 
and where the right- and left-invariance are enforced by sandwiching the product 



<s 



between two Haar intertwiners. This Haar projector is just the group averaging in 
(|17p translated into the spin representation. The factorizing property means that P e 
can be written in the following way 



P. 



PHaar ' P e ' PHaar = PHaar ' [ E pi ® E p2 



' ' ' Ep n ) ' PHaa 



(18) 



which again is (jTT|) group-Fourier transformed (see figure [3]). Here E Pi are (usually 
not covariant under the group action) maps acting on pi and with matrix components 
given by the group Fourier transform of the functions E, i.e. 



(E p r n = £ e p kdd ,I(p,kr t (lHp,k))\ 
using the notation from (fTU|) . 



(19) 



:.d',i 



p 



— E — — 



E 



Figure 3: Graphical notation of the relation between the OSFM version of the edge 
propagator P e and the one used in ©, which is also used in the BC and 
EPRL models. 



4.2. Group Field Theory formulation 

The Group Field Theory (GFT) formulation [33l [34l EH [36] of the spin foam state sum 
is based on a quantum theory of fields defined on the group G, such that the Feynman 
graphs are Spin Foam amplitudes associated to 2-complexefl For one of these 2- 
complexes the group holonomies are pairs of holonomies g e v f-,g%p where v C e, e' C /. 
Both are interpreted as going from the vertex v to the center of the face /, and each is 
associated with one of the two edges e, e' that border / and meet at v (see figure 0]). 
As usual, we denote the inverse by (g^j)^ 1 =■ g%- A natural class of models arises 

if one indeed imposes g^ = g*j, as Figure 3] suggests, but the GFT models allow for 
relaxation of this condition. 

The ingredients for the GFT amplitude associated to C are edge propagator functions 
II e , each depending on all g^ v for a given e, and vertex functions V v , each depending 



1 Where the specific form of the action is designed to control the class of appearing 2-complexes, see 
e.g. [37]. 
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Figure 4: Variables in the group field theory formulation. 



on all for a given v. The functions have to satisfy the following symmetries: 

Vv(---,gZf,9vf,---) = V v {...,ggl f ,ggi } ,...) (20) 

n e (-.-,^/,fl^/,..-) = n e (..., g% f g,g*, f g,...) 

The symmetries (|20j) make sure that V effectively only depends on the combinations 
gf v g%f, and the propagator functions IT depend only on g^gj v ,. Also, both functions 
V v and n e have to be gauge-invariant with respect to gauging at an edge, i.e. g e v ^ — > 
<7 for fixed v, e, and all / which appear as variables in the respective function. As 
a result, the II e will be invariant under left or right shift of its arguments. The V v 
can be decomposed into spin network functions of the boundary graph of the vertex 
V, which gives the usual spin representation of the vertex amplitudes. 
The state sum for the GFT models is given by 

z gft = [ dg% s If II, (>/;., g%, Sn ) [] v - (C h » • • • » 9%0 ( 21 ) 

•* e v 

The i?-function form ((5]) for the state sum Z can be written as a GFT state sum (|2"Tj) 
for u = 6. If the product of the h e f around the faces are constrained to be the unit 
element, one can write (uniquely up to global right shift of all g%f for some fixed / by 
the same g 6 G) 

he! = 9 e vf 9%', (22) 
imposing g^ = g^ in the state sum. This leads to 

Z = I d9 e vf ]JC(g e vf g e fv „...)l[5(g e vf g}' v ) (23) 

= J d 9i f n M^-X'fJUMglh-^J 



with vertex function 



M9^---'9 e v ?J ■= II 'H.'/',/.'//, ) (24) 

fDv 
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and propagator 

n e (gl h ,...,gl, fn ) := C(g e vf g%„...) . (25) 

4.3. Vertex amplitude formulation 

A formulation closely related to the GFT formalism can be given in terms of holonomies 
labelled by wedges, i.e. by a face /, and pairs e, e! that meet at a common vertex (see 
Figurc[5|). The holonomy can be thought of as going from the center of e' to the vertex 
v, and from there to the center of e. 




Figure 5: Holonomy variables in the vertex amplitude formulation. 



One can arrive at a formalism using the g[ by imposing g e > v g ve = g[ (as suggested by 
figures [TJ and [5]) in the state sum and integrating out the gauge holonomies g ev , 
arriving at 

^vertex = / dg ei dg{ J[ U (Gf) [] E( 9ef ) J[ W v , . . . , gfr) (26) 

f ef v 

where Gf = glg e fg{i • • ■ 9e'f is the ordered product of group elements around a face 
/, according to its orientation. The vertex function W v for a vertex v is given by 

WM 1 , • • • , 9t) = / dg ev I] S{g efv g w { 9 i)- 1 ) . (27) 

fDv 

The vertex functions W v have a symmetry that corresponds to gauge transformation 
associated to edges, i.e. for an edge e beginning at vertex v, V v is invariant under the 
transformation of its arguments 

f glh e if f D e 
d ■+ t 7 , (28) 

This invariance allows to interpret W v as a G-spin network function for the boundary 
graph of the neighbouring 2-complex of the vertex v. A very similar formulation that 
uses G'-spin networks instead was given in [23] . 
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Writing the state sum model in the form (|26|) allows for generalizations of the model, 
by relaxing the condition (|27[) . allowing for more complicated interactions at the ver- 
tex. These non-trivial vertex functions actually do appear in a renormalization group 
flow via the n — 1 Pachner moves of state sum models, which is why using the space 
of all models of the form (|26p allows for computing an exact renormalization group 
flow in this case. The cost of this is that this space is much larger than the space of 
E- function models, using state sums of the form ((6]). Still, the example in section [6. 3. II 
shows that truncating the flow to this subspace can still capture essential features of 
the phase space, as well as its fixed point structure. 

The form (|26|) can be simplified by absorbing the iJ-functions into the vertex func- 
tions W v . To this end, we write the E- functions as convolution of two functions, 
i.e. 



E(g) = J dkF{k)F(gk). (29) 

The function F has the same symmetries than E, i.e. F(g) = F(g' g(g')^ 1 ) = F(gh) 
for g,g' € g, h 6 G'H Thus we get one more integration over G for each pair ef, 
resulting in 




Figure 6: Wedge holonomies in the alternative vertex amplitude description. 



Performing the variable transformation g e f —> g e fk e f and g[ — > h( = k e fglg e ,^ and 
carrying out the integration over the g e f and k e f, one arrives at the form 



^Vertex 



J dhl\[u{G f )X[V v (hl...) (31) 



2 It should be noted that _F is not unique, but depends on choosing a collection of phase functions. 
If, however, E is positive, then F can be chosen positive as well. 
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where G/ is the ordered product of the k[, for all v C /, according to the orientation 
of /. The new vertex functions V v are given by 




The vertex functions V v inherit properties from both W v and F. Firstly, it can be 
readily seen that V v can - just as W v - be interpreted as a G-spin network function 
on the boundary graph around the vertex v. Furthermore, the right-invariance of F 
under the subgroup G' implies that V v is invariant under the symmetry 

hi -+ hih f v h 2 , h!,h 2 G G'. (33) 

Hence V v is not only a function on elements of G, but rather on G = G'\G/G'. If 
G' happens to be a normal subgroup of G, G in itself is again a group that carries 
a right and left G action. In this case V v can be interpreted as function on a graph, 
depending on G-elements on the edges, invariant under gauge transformations at the 
vertices with elements from G. 

It should be noted that one can arrive at the same version of the state sum (|3ip by 
starting from a version of G e -functions and vertex functions W v . Since the G e -functions 
are more general than the -©-functions, in this case the resulting vertex functions V v are 
just G-spin network functions, without any specific properties regarding the subgroup 
G'. It is this form, i.e. equation (|3ip with general G-spin network functions V v and 
class functions u>, which is the most general of all the presented formalisms. 



State sum 


symbol 


variables 


parameterizing functions 


E- function 


Z 


9ev, 9ef 


LJ, E 


OSFM 


-^OSFM 


h ef 


w, G e 


GFT 


Zqft 


9%f 


n e , V v 


Vertex function 


-^Vertex 


gt, 9ef 


u, E, W v 


Vertex function (alternative) 


-^Vertex 


H 


w, V v 



Table 1: Different formulations of the Spin Foam state sum, using different variables 
and parameterizing functions. 



5. The BF, BC and EPRL models 

Here we will shortly describe the structure of the current spin foam models [51 [T51 [TO] 
in the holonomy language. To do so we specify now to G = Spin(4) and G' = SU(2) = 
Spin(4)diag the diagonal SU(2) subgroup, with irreps labelled by p and k respectively. 
Also, the face weights arc usually set to be the delta function on the group. 

The operators E described in section 14.11 for the various models are given by the 
following, see for instance [2]. 
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• For BF theory we have 

E BF = 1, (34) 
which reduces the edge projector P e to be the Haar projector. 

• In the Barrett Crane model the maps E p which act on representation space p, 
restrict to representations p = (fe, k). Furthermore E p just involves the injection 
maps (described in section l3~Tj) I(p, 0) from p to the representation k = 0: 

E* c = d e (p)J26(p,(k,k))I(p,0)I(p,0y . (35) 
fe 

where d e (p) is an unspecified edge measure factor. 

• For the EPRL model the maps E p also restrict to a fe-parameter family of rep- 
resentations, which is now of the form p 1 {k) = f^p-A, - ^ k^j with 7 being the 
Barbcro Immirzi parameter: 

E^ PRL =d e (p)^26(p,p,(k))I(p,k)I(p,k^ . (36) 
fe 

The injection map picks out a representation k p which is determined by p via 
the relation p = p 1 (k p ). 

• The FK model is constructed using coherent states [35]: a description of the 
corresponding edge operators and functions can be found in [2j [28] . 

The edge functions E(h) can be obtained from the edge operators E by Fourier 
transform and using the parameterization in section f3. II 

• For the BF model we just have E BF (h) = 80(h), i.e. the E function is given by 
the delta function on the group. 

• For the Barrett Crane model we obtain 

E BC (h) = MP) dim(p) S (p, (k, fe)) tr, (D p (h) I(p, 0)I(p, 0)t) . (37) 
p.k 

For the choice d e (p) = 1 this results into 

E BC (h) = Ssuwih+h-- 1 ) = e Spin{4)diag (h), (38) 

where, for a general subgroup G C G we write 

Q G '(h) = [ dtiSihh'- 1 ) (39) 
Jg> 

for the delta function that forces a group clement to lie in the subgroup. Note 
that the i?-function has an enhanced symmetry with respect to the subgroup 
G' = Spin(4)dj aff , it is not only invariant under conjugation with elements in G' 
but also invariant under left and right multiplications. This explains some of the 
special properties of the BC model as compared to more general models. 
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• For the EPRL model the function E EPRL (h) is given by 

E EPRL (h) = dim(p) S(p,p^(k)) tr p (D p (h) I(p, k)I(p, fe)t) . (40) 

/.'.A- 

As for the BC model this will be a distribution, with singular support described 
below. 

The E functions arc distributional for all the cases described here and have support 
on the following critical manifolds: 

• BF: 1 

• BC: S\J(2) diag 

• EPRL: M 7 = {g G Spin{A)\3L G su(2) : g = (cxp(i), cxp( 7 L))} . 

For the EPRL model this is shown in reference |29j . For our purposes it will be 
convenient to characterise the EPRL support in terms of the group. This can be done 
for example by choosing a U(l) subgroup of U(l) 7 <E Spin(4) generated by (L,-fL). 
The critical manifold is then simply S\J(2)di a g >ad U(l) 7 , where t> a d is the adjoint 
action. 

5.1. General (finite) groups 

We can now mimick the three types of spin foam models defined above for finite 
groups. To do so we assume that the most important part of the amplitude is the 
divergent, distributional part. To define models for finite groups we choose functions 
with support on the appropriate subsets of M C G, that is by replacing the delta 
functions by normalized 0m (<?) functions that arc defined to be pg^Vg G M and 
otherwise. The _E-functions are then: 

• BF: E(h) = e (/i) with e being the trivial subgroup consisting of the identity 
element. 

• BC: E(h) = Q G ,(h) 

• EPRL: Ec(h) = <dm c (h), where C is a cyclic subgroup of G and Mq = G'\> a dC. 
This generalizes the BC and the EPRL model to arbitrary groups G' C G. 

5.1.1. Example: S3 

Let us consider as an example the finite group G — S3, the group of permutations of 
three elements, as gauge group. It can be written as 

S 3 = (a,b\a 2 = b 3 = {abf = 1). (41) 

The group has six elements {e, a, 6a6 _1 , b 2 ab~ 2 , b, b 2 }, where e is the unit element, 
{a, 6a6 _1 , b 2 ab~ 2 } are odd permutations (two-cycles) and {b, b 2 } are three-cycles. 
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The group has subgroups isomorphic to Z2, generated by either a, ab and ab 2 . 
Furthermore there is a Z3 subgroup e,b,b 2 which is normal in 5*3. To specify the 
models further we have to choose one of the subgroups and consider E functions 
which are invariant under the adjoint action of the subgroup G on G. The orbits 
under the different subgroups, including the group G = S3 are 

G' = S3 : {e}, {a, bah" 1 , b 2 ab~ 2 }, {b, b 2 }; 
G' = {e, a} = Z 2 : {e}, {a}, {bab" 1 , b 2 ab~ 2 }, {b, b 2 }; 
G' = {e, b, b 2 } = I3 : {e}, {a, bab~\ b 2 ab~ 2 }, {b}, {b 2 } . (42) 

Choosing G' = Z2 as subgroup, we have a priori a four parameter space of E- 
functions 

E(g) = 13,6(9, e) + p 2 S(g, a) + p 3 (S(g, bab' 1 ) + d(g, b 2 ab~ 2 )) + (3 4 (6(g 7 b) + 6(g, b 2 )). 

(43) 

Taking some normalization condition into account, we obtain three parameters. As 
will be explained in section I6.lt models in which the E function is invariant under 
conjugation of the full group can be rewritten as standard lattice gauge theory models, 
i.e. the E functions can be absorbed into the face weights u>. This will be the case for 
Pi = 

The BC model for the choice of this subgroup Z2 corresponds to /3x = Pi = 1 and 
P3 = /?4 = 0. Hence it will in general define a non-trivial model, which cannot be 
rewritten into a lattice gauge theory. 

We can define an EPRL like model with Z3 as cyclic subgroup, which leads to 
f3i = /?4 = 1 and /?2 = /?3 = 0. This model is also a Z3 subgroup model (the BC model 
with subgroup Z3) and actually rewritable as a standard lattice gauge theory. For this 
choice of parameters we have 

E(g) = e Z3 (g) (44) 

and as Z3 is a normal subgroup the function E will be invariant under conjugation. 
Indeed the BC model with a normal subgroup will define a topological theory which is 
of BF type on the quotient group (here Z 2 ) and of high temperature type (i.e. E is a 
constant function) on the subgroup itself. We can define two different one-parameter 
families which connect this topological theory to either the high temperature/strong 
coupling fixed point of the full group (Pi = P2 = P% = Pi) 

E(g)=Q Z3 (g)+jQ odd (g) (45) 

(here o dd = 1 on odd permutations and vanishes on even permutations) or to the BF 
theory on the full group 

E(g) = 6(g,e) + p 4 (6(g,b)+6(g,b 2 )) . (46) 

Another example to define a BC model more closely connected to the SU (2) x SU(2) 
model is to take G ~ S3 x S3 as gauge group and the diagonal group as subgroup. We 
will discuss this model shortly in section 16.31 
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6. Coarse graining 



In this section we will consider coarse graining of the holonomy spin foam models 
in dimensions two and three. To this end we will utilize in the first part techniques 
from standard lattice gauge theory, that is use the similarity of holonomy spin foams 
to standard gauge theories. Indeed, as for gauge theories, coarse graining in two 
dimensions can be done exactly. The holonomy spin foam models feature in addition 
to the face weights, which are the central dynamical entitles for standard lattice gauge 
theory, the edge weights or edge functions E, which will lead to additional factors in 
the coarse graining formulas. 

In three dimensions certain models, including BC type models, can be rewritten into 
models of standard lattice gauge theory type. We will identify this class of models and 
coarse grain these models by applying the Migdal Kadanoff scheme. 

These results can be compared to another type of coarse graining scheme, which 
involves a hierarchical lattice on which coarse graining can be performed by applying 
repeatedly 4—1 Pachncr moves. This allows to do an exact coarse graining. In 
particular we will test whether the paramctrization of the holonomy models with E 
function captures the relevant dynamical data. That is although the general coarse 
graining flow leaves the model space described by the E functions we can project back 
to this space and test whether we still find the fixed point structure of the exact flow. 

Attempting concrete (numerical) calculations we face the difficulty that the models 
are not necessarily finite for (compact) Lie groups. To avoid these difficulties numerical 
investigations will be (mostly) performed with finite groups, in particular S3, for which 
we defined a class of models in the last section. 

6.1. Coarse graining of edges 

Let us consider the formulation of the models (fT4"|) . where to every edge we associate a 
C-function and to every face face weights uj. One important difference of these models 
to standard lattice gauge theories is, that in general, these are not invariant under edge 
subdivision, or the inverse operation, coarse graining of edges by removing two-valent 
vertices. 

Therefore wc will briefly discuss the coarse graining of edges. This is most conve- 
niently done in the spin basis where the C-function is expressed as follows 

C(hi,h 2 , - ■ - h n ) = 2J dim pi Cp lt ... tPni dd< (i>d/ | D Pl (fti) ® • • • <g> D Pn (h n )L d ) . 



Here {i^} is an orthonormal basis of the intertwiners and (• | •) is the inner product on 
L 2 (G"). 

The coefficients C can be computed from the coefficients e k pdd , defined in (fl"0|) if the 
C-function is derived from an _E-function as in 

C(h efl , . . . , h efn ) = f dg L dg R J] dh efi Y[ E(h efi ) Y[ 5(g L h efi g R h- f ) ) . (48) 
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Expanding the 6 function 5(g) = dim px P {g) and using the identity 

/ dg L dg R J[ Xp , (9l h e f z 9rKu ) = 
J i 

Y, tr Mi* i ® d p< (^/« ) i ^ ^ i ® D P' (Kb ) ( 49 ) 



we obtain 



C pl ..., Pndd , =Hdh efi ]jE(h efi ) U d \(g)D Pi (h- f \)i d , 



(50) 



Using the explicit form for the E function (fTU|) yields 
C 



..., Pn dd'= Y[^. d . d ,_(i d \^}I(puh) di I^(p i ,ki) d! . 

ki,did' i i 



hi' 



(51) 



With the form (|T71) of our edge function C at hand we can now consider the coarse 
graining, where a subdivided edge, i.e. two edges e, e' joint by a two-valent vertex, is 
coarse grained into one edge n. All edges bound the same faces /i, . . . , f m , see Figure 
UJ As C-functions are composed by delta-functions, the C-function for the new edge 
n is given by 



c-w {h nfl , . . . , h nf j = J n dA eA n d/i -'/. 



<?W (h efl h efm ) (h e , fl , . . . , heffj J] A ) 



(52) 



Expanding all functions involved as in (|47p and using the orthogonality of the basis i d 
we arrive at 



C 



(n) 



pi,--- ,p n ,dd" pi ■ ■ ,p n ,d"d' 



(53) 



This gives the renormalization flow equation for the coarse graining of edges by re- 
moving two-valent vertices. Note that the flow equations for this operation are local. 




Figure 7: Removing a trivial subdivision of an edge. 
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In 2D and many 3D models the C-function simplifies drastically, so that the models 
can be reduced to standard gauge theory. To explain how this happens, we define an 
n-valent C-function to be factorizable, if it can be written as 

C(h efl ,h ef2 , . . . ,h e f„) = J dg L e(g L h efl ) e{g L h ef . 2 ) ■ ■ ■ e(g L h ef J . (54) 



where e is a class function. Let us notice that every function of the form (|54p can 
be rewritten into a left and right group averaging of i?-functions since due to the 
symmetries of e we have 

C(h efll h eh ,...,h efn ) = J &g L dg R e(g L h efl g R ) e(g L h eh g R ) ■ ■ ■e(g L h efn g R ) , (55) 

hence E = e. The form introduced above allows us to integrate out h e f and to absorb 
the e factors into the face weights w. If we define the effective face weights 

u)' f = J Y[ dh e fe(h e f) uj f {- ■ -gvehefgeyi ■ ■ ■ ) , (56) 

eCf 

these are again class functions of the group element • • ■ g ve g e v' • • • ■ (Note that the face 
weights ui'j might depend on the face /, in particular if faces differ in the number of 
adjacent edges.) Explicitly, using the Fourier transform 

cy(. 9 )=]TdinW P Xp(3) (57) 

p 

e{g) = Y,&™Pt P X P {9) (58) 
p 

we obtain 

QP = ( £ P)«(eC/) ~p_ (5g) 

In this way every spin foam with factorizable C-functions can be rewritten as a stan- 
dard lattice gauge theory. This factorization property holds for every two-valent edge 
and in the case of BC/EPRL/FK also for every three-valent edge. 

If we do not integrate out the h e f and keep the C-functions, the Fourier transform 
of these will be given by 

c Pu ..., Pn , dd > = (ll £p '^ Sdd ' ■ ( 6 °) 

This will simplify the renormalization flow (|53[) . 

6.2. Coarse graining in two dimensions 

Here we will consider a 2-complcx for which all edges are 2-valent, i.e. bound exactly 
two faces. In this case the coarse graining flow can be computed exactly. The procedure 
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and result is quite similar to real space renormalization of 2D lattice gauge theories 
[39] . As pointed out in section 16.11 compared to standard lattice gauge theory one 
obtains however an additional renormalization flow from the edge weights. Starting 
with the form of the partition function 



Z = J fn^'l (n^eMe/oj (j["f( h ef h e'f---)\ > ( 61 ) 

we have to discuss the renormalization flow of the C-functions and the face weights w/ 
The C-function is by construction invariant under left and right multiplication (in 
both entries). Hence it can be expanded into a basis of gauge invariant functions, that 
is 

C(h u h 2 )^Y, d P d[m (PM h ^ 1 ) ■ (62) 
p 

Such a C-function is factorable in the sense of the Fourier expansion of the 

factorization property gives 

C(h!,h 2 ) = [ Age{gh 1 )e{gh2) = ^ dim^Xp^i^ 1 ) . (63) 
J p 

Hence we just need to choose the Fourier coefficients e p such that e p e? = C p is satisfied. 
The coefficients C p can be expressed using the parametrisations through e p k dd , as 

k H dd 1 

where we used that e? = e^ d , d , which can be shown to hold in general |28j . That is 



"kdd' 

from the choices e kdd , for the basis coefficients of E only the contraction matters 
in a 2-dimcnsional theory. 

The renormalization flow equations for the coarse graining of edges (|53p are in this 
case given by 

£(n) = cV>c&) . (65) 

The other coarse graining step that will appear in two dimensions is integrating 
out the variables over one edge e that divides two faces /, /', see Figure [5] The face 
holonomics starting at the source of e will be denoted by p/ = h e fHf and <?/' = 
hefHfi. The new face weight is given by 

An) 



f iE f H?) = j\u,, j „,, r . j U,, J U / ^ r U,, r U / , ( V, f .>,^ , (66) 
which for the Fourier coefficients gives 

G$f = a,{p)Q r (p)C^ . (67) 
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In two dimensions all coarse graining steps can be decomposed into the two oper- 
ations we discussed, that is the inverse of subdividing a face with one edge and the 
inverse of subdividing an edge with one vertex. This allows us to write the complete 
partition function of a 2-dimcnsional complex as 

Z = J2("(P)) F (C P ) E (68) 
p 

(assuming that the face weights and edge functions are initially the same for all faces 
and edges). 

The recursion relations (|55|) and (|57)) enable us to discuss completely the renormal- 
ization flow in two dimensions. Starting with a 'heat kernel' (C p is the Casimir of 
P) 

C p = cxp(-/?- 1 C p ) (69) 

for the edge weights, we flow from zero temperature j3 = oo to infinite temperature 
(3 = 0. Note that zero temperature j3 = oo, C p = 1 corresponds to a C-function 
C(hi,h 2 ) = S^ih^ 1 )- This converts the partition function into one for standard 
lattice gauge theory. Choosing non-zero temperature corresponds to a regularization 
(introducing a non-trivial propagator for the edges) for instance in group field theory 
[35] , One should however be aware that this choice leads to a flow into the high 
temperature fixed point, both for the edge factors and the face weights. In this fixed 
point only the values for the trivial representation contribute (assuming C p= o = 0). 

The renormalization flow for the face weights behaves similarly, with the addition 
that it can be set off from the zero temperature fixed point u>(g) = 5(g) by a non-trivial 
choice for the edge factors alone. Note that the zero and high temperature fixed points 
are not the only fixed points, as will be obvious from the discussion below. 

In higher dimensions the subdivision of an edge by adding a two-valent vertex, 
as well as the subdivision of a face by inserting a two-valent edge can be argued to 
be operations under which the partition function should remain invariant [40j . Indeed 
this notion has been connected to the anomaly-freeness of the path integral measure in 
spin foams with respect to diffeomorphisms [4lJ [TUJ 02] and has the advantage of fixing 
uniquely the edge and face measure factors [40j . Such an invariancc would require C p 
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and oj(p) to take values either equal to one or to zero. Any choice would lead to a fixed 
point in the two-dimensional theory. Also note that for the two-dimensional theory 
it is equivalent to either choose uj(p) = or C p = for a given p. That is all C p can 
be chosen to be equal to one. Hence we just need to consider the lattice gauge theory 
parameter space, reduced to the choice oj(p) = 9(p) where 8 takes only values and 1. 



6.3. Dual simplicial lattices in three dimensions 

We have discussed that the renormalization flow for the two-dimensional theory agrees 
with the one in lattice gauge theory in case that one adopts the condition of invariance 
under edge subdivisions. This conditions is also a reasonable requirement to hold for 
higher dimensions and higher-valcnt edges. Indeed the Barrett-Crane model can be 
easily made to satisfy this condition. For the EPRL model the construction [121 [35] is 
more complicated and leads to a model that falls not anymore into the description by 
i?-functions (would however be still in the more general class of models parametrized 
by C- functions). 

Let us assume that we can construct a spin foam model where such an invariance 
holds. How much freedom is left for the parameters in this spin foam model? In 
a dual simplicial lattice in three dimensions the edges (dual to triangles) would be 
three — valent, i.e. adjacent to always three faces (as triangle has three edges). The 
C-function associated to such three-valent edges is a function of three arguments, 
invariant under left and right group multiplication. Hence it can be expanded into a 
basis of gauge invariant spin networks, which are based on three edges with common 
source and target vertex, see Figure O 




Figure 9: The structure of the C function for an edge shared by three faces. 



These spin networks are labelled with representations at the edges and with inter- 
twiners between the three representations, that is we have 

C(/ii,/i 2 ,/i 3 ) = 2J ~~[ dim p t C Pl 

,P2,P3,dd' (I'd') D Pl (hi) ® D Pl (h 2 ) ® D Pn (h 3 )od) 

pi,p2,P3,dd' i 

where labels orthonormal basis of invariants. We will consider groups, such as 
SU(2), SO(A) and the finite group S3 for which the intcrtwincr between three repre- 
sentations is canonical, or in other words for which the tensor product of two rep- 
resentations is multiplicity free. In this case the intertwiner labels d, d! in (|70[) can 
be omitted. Also the sum in (|70l) is only over triples of representations admitting an 
intcrtwincr. 
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Gluing two edges with a two-valent vertex leads to new C-coefficients 

Mn) = cM') r<( e ') (H\\ 

Pl,P2,P3 Pl,P2,P3 Pl,P2,P3 ^ ' 

Demanding invariance under edge subdivision fixes the coefficients to be equal to 
one or zero, depending on the three representation labels. In the BC and EPRL model 
(with G isomorphic to G' x G' group) this condition even factorizes as in equation 
([50]). that is we have 

Cp 1 ,p 2 ,P3 = e Pl e P2 e p3 (72) 

where e p takes only the values zero or one. One might expect such a factorization 
property to hold from the construction of the edge propagator C as group average 
of a product of iJ-functions. Note however, that even in this class of theories the 
factorization property docs not hold in general. For a counter example see appendix 
IA1 Also the BC model for 5*3 with Z2 as a subgroup provides a counter example. 
As one can check explicitly, the corresponding three- valent edge function C(h\, /12, ^3) 
does not lie in the subspace generated by the conjugation invariant E'-functions, which 
define factorable models. 

The following considerations restrict therefore to the class of theories where the 
coefficients factorize as in (|72|). In this case we can absorb the choice function e p 
into the face weights w / , whereas the C can be chosen to be the same as for standard 
lattice gauge theory (where invariance under edge subdivision trivially holds). In other 
words we are reduced again to consider the class of theories described by standard 
lattice gauge theories. This motivates us to test approximation methods developed 
for lattice gauge theory, notable the Migdal-Kadanoff approximation [SSJJ [H] . There 
is one draw-back which is that the Migdal-Kadanoff approximation was developed 
for regular (cubical) two-complexes and for hierarchical lattices |45j (where it can be 
made to be exact). That is we have to assume that the difference of a simplicial versus 
a cubical lattice does not influence the large scale dynamics. 

The Migdal-Kadanoff approximation is a truncation to local couplings for the real 
space renormalization flow of lattice gauge theories. As for the flow in two dimensions 
it results into a recursion relation for the face weights w/. In the version of |39j we 
have for the coarse graining of a three-dimensional cubical lattice under the doubling 
of the lattice constant 

~(n) \4,~ ,4 dim (Pl) dim (P2) p ,„ Q x 

w p - d~ 77) ■ ( 73 ) 

P1-P2 " ' 

Here n p pi p2 is the number of copies of representations p in the tensor product of p\ 
with pi. For the case we are considering here, this number is either one or zero. 

Note the essential difference of ([73]) to the two-dimensional recursion relation 
which is that the coefficients of different representation labels are coupled with each 
other. In the BC and EPRL spin foam models for a group G = G' x G' only a 
(one-parameter) subset of all possible representations is allowed initially, however the 
recursion relation (|73p will spread the range of allowed values. 
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The implementation of (|73|) for a finite group is straightforward. Here we consider 
the group S3 x S3 with a BC like amplitude. There are three representations for 
S3 which we denote by j = 0,1,2. The representation j = 2 is two-dimensional 
whereas the others are one-dimensional, given by the trivial and the (permutation) 
sign representation. For S3 x S3 we have representations (j + , j~) with j + , j~ = 0, 1, 2. 
The initial face weights are given by 

<Vj- = S 3+,j- > ( 74 ) 

which we put into the recursion relation ([75]) . The coupling coefficients are given by 

nf+'Q r+ .-,=n j+ + . + n r - ._ (75) 

where ■ is invariant under permutations of the three labels and equal to one for 
the following combinations of representation labels 

(0,0,0), (0,1,1), (0,2,2), (1,2,2), (2,2,2) . (76) 

For all other combinations the coefficients are vanishing. After each renormalization 
step we normalize the amplitude so that <I>o,o = 1. The result is that after ca. 10 
iterations the theory is at a fixed point 

Qj+j- = (S j+fi + Sj+,i){8j- j0 + Sj-,i) , (77) 

These weights correspond (modulo normalization) to a face weight 

^(91,92) = @even(gi)@even{g2) (78) 

which defines a BF(^2) theory, where this Z2 can be understood to be the the quotient 
group of the diagonal subgroup with respect to the Z3 subgroup of even diagonal 
elements. 

Let us also discuss the SU(2) x SU(2) BC model. Representations are given by 
(j + where the spins can take integer or half integer non- negative values and the 
dimension of the representation is given by (2j + + l)(2j~ + l). The coefficients 

are - a is equal to one, if the triangle inequalities ji + jj > jk arc satisfied and if the 
sum ji + J2 + J3 is an entire number, otherwise it is zero. The initial face weights are 
again given by (f71| . 

However the sum (|73[) diverges. Therefore we introduce a cut-off and sum values 
only up to j < A. Again we normalize after each iteration step, such that 0)0,0 = 1- 
An impression of the behaviour of the flow can be obtained from figure 16.31 

As one can see the face weights flow very fast into the high temperature fixed point 
and the result does almost not vary with the cut-off. Note also that whereas initially 
only a one-parameter subset of face weights is non-zero, this range is spread during the 
recursion process. Therefore one can also expect a similar flow for the EPRL model, 
where the initial face weights are also only supported on a one-parameter set of labels. 

Of course the introduction of the cut-off can be questioned. On the other hand 
this result is supported by a general theorem, stating that in 3D gauge theories with 
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Figure 10: This figure shows the face weights u>j+ j- as a function of 2j + , 2j~, which 
label the x and y axis respectively. The initial configuration Qj+j- = 5j+ j- 
of face weights is shown in blue. The red dots give the configuration after 
one iteration, where we show only a dot for values of uij+ j- > 10~ 4 . The 
cut-off is A < 15. As can be seen from the picture the initial diagonal con- 
figuration is spread to non-diagonal arguments, but the values for the face 
weights are decreasing fast. After two iterations (yellow- green dots) a large 
portion of the ujj+j- is already smaller than 10 and the configuration is 
almost at the high temperature fixed point. 

groups U(l), SU(N) flow under the Migdal Kadanoff relations to the strong coupling 
/ high temperature fixed point gSJ 05] . 

Now, we made several approximations, which could be questioned. The first is to 
replace a simplicial lattice with a regular one, the other is the Migdal-Kadanoff approx- 
imation itself. To support the findings here, we considered an exact renormalization 
flow for a special simplicial lattice, which will be discussed more extensively in the 
next section. This lattice is obtained by iteratively applying 1 — 4 Pachner moves to 
a tetrahedron, see figure [TT] 

Such a 1 — 4 Pachner move subdivides one tetrahedron into four. To these four one 
can again apply the Pachner moves obtaining sixteen, and so on. Going backwards, 
we can coarse grain blocks of four tetrahedra into one. The advantage of such a 
hierarchical lattice is that the renormalization flow stays in the set of local couplings. 
Truncations, which are usually unavoidable in real space renormalization due to the 
generations of non-local couplings, are not necessary. This renormalization flow on 
the simplicial lattice is most convenient to consider in the spin representation with 
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Figure 11: The 1-4 Pachncr move can be iterated to generate an actual flow in the 



vertex weights. Here one is actually working with the space of vertex weights for an 
S3 model. In the spin representation the BC model with S3 x S3 is included in this 
space. The flow of this model terminates in a hxed point describing a BF(Z2) theory, 
which confirms our findings with the Migdal Kadanoff truncation. 

Another example, which we can consider within G = S3 is with weights u>o = 1 
and £j2 = lH With the Migdal Kadanoff recursion relation (|73|l it flows to BF theory 
for the group S3. With the hierarchical lattice flow, explained in the next section we 
obtain a different result: here it flows to the strong coupling/high temperature fixed 
point corresponding to have only the coefficient Qq = 1 to be non- vanishing. 

This difference can be expected already due to the different structure of the Migdal 
Kadanoff flow and the hierarchical flow. A typical feature of the Migdal Kadanoff flow 
is that the u> are convoluted with each other, which leads to a spreading effect of a) as a 
function of the representation labels. For the S3 example it is therefore possible to flow 
from a configuration with Q\ = 0, i.e. where the sign representation does not appear, 
to the BF fixed point, in which Cb\ = 1. In the flow defined by n — 1 Pachncr moves 
this does not happen. The flow is described by vertex functions which depend on (six) 
representation labels. If initially the vertex functions are zero for any of the labels 
giving the sign (or any other) representation, this will be preserved by the flow. This 
is however a special feature of the n — 1 Pachner moves. To get a more complete pic- 
ture of the flow for more general lattice one has also to consider other Pachner moves, 
which will require some implementation of a truncation. The Migdal Kadanoff method 

3 This model corresponds to the choice of parameters @i = 5, /?2 = P'i = lj/84 = —1 in d 43 1 ) and 
ei = 0, e20 = e2l = 1 in the next section. 



space of models. 



26 



is one particular truncation available for lattice gauge theories. Other methods, such 
as tensor network rcnormalization |47l I15j are more general and in principle also appli- 
cable to spin foam models [TS]. Tensor network renormalization can be motivated by 
using Pachner moves |47j , and indeed the Pachner moves different from the n — 1 move 
are the ones where the truncation is implemented in tensor network renormalization. 
This choice of a truncation can be also encoded into embedding maps which are closely 
related to the concept of cylindrical consistency in loop quantum gravity |48) . 

6.3.1. Renormalization group flow of an S3 model via 1 — 4 move 

In the following we will consider the flow of 3D models with group S3 generated 
under the 1 — 4 move. We described this group in section 15.1.11 We also choose 
G' = {l,a} — Z 2 as subgroup. 

We will define the model on the 2-complex ki dual to a tetrahedron. Its topology 
allows to glue four of them together to form the 2-complex K2 dual to four tetrahedra 
which are the result of a 1-4 move. The boundary of that 2-complex will be the same as 
the one dual to one tetrahedron, apart from trivial subdivisions of edges. The effective 
model on this, arising from integrating over internal degrees of freedom, can then be 
compared to the original model, by evaluating Z Kl and Z K2 on the same S^-boundary 
spin network functions. Iterating this process generates a flow in the space of models. 

If we define the model in terms of vertex functions, the flow can be carried out 
exactly. This is a special feature of all l-(n — 1) moves, since they do not generate 
any non-localities. These moves have another interesting property, which is present 
for all models and dimensions: The vertex function can be interpreted as a function 
depending on spins and intertwiners. If one starts the flow with a function that is 
constrained to vanish whenever a specific spin or specific intertwiner is present, then 
this constraint is conserved under the l-(n — 1) move. Therefore, if certain degrees 
of freedom are switched off, they stay switched off during the flow. This observation 
allows for a very convenient investigation of the RG flow, since it is automatically con- 
strained to lie in certain submanifolds of the model space, all of which intersect in the 
high temperature fixed point. Note that this is a direct result of working on a hierar- 
chic lattice, because the flow does not generate any non-localities. It is quite different 
from the hypercubic lattice flow described in the last section, where non-localities are 
truncated by the Migdal-Kadanoff approximation. 

Let us turn to £3 once again. As previously described there are two 1-dimensional 
representations (denoted 0s 3 and ls 3 ), and one 2-dimensional representation, denoted 
by 2s 3 . Since 0s 3 is the trivial representation, the only nontrivial tensor products are 
given by 

1 S3 ®ls 3 = S3 , l Sa ®2s 3 = 2 Sa 

2 S3 ® 2 S3 = S3 © Is, © 2 S3 

The group Z2 has two 1-dimensional representations, denoted by 0z 2 and lg 2 , satis- 
fying lz 2 (§3 lz 2 = 0^ 2 . Decomposing the representation spaces of 5*3 with respect to 
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TL-i yields 



0s 3 = Z2 

Is, = lz a (79) 
2s 3 = Z2 © 1 Z2 

After gauge-invariance and discrete symmetries of the tetrahedron, there are eleven 
gauge-invariant combinations of 53-irrcps along the edges of a tetrahedron. In other 
words, a vertex function is given by specifying the eleven different values of the 
symbol for 53 

If we specify a model in terms of i?-functions, then there are fewer parameters: Due 
to the decomposition (|79p the most general _E-function is - after normalization - of 
the form 

E{h) = 1 + eiBgn(fc) +2e2oPa s ,Woo + 2 e 21 p 2s3 (h) n (80) 

where sga(h) = trp ls3 (h) is the signum of the permutation h € S3, and (00) and (11) 
specify the upper left and the lower right entry of the 2 x 2-matrix p 2s (h). A model 
is then given by the three parameters ei, e2o and e2i0 

Note that every model in terms of E'-functions can be rewritten in terms of vertex 
functions, but since these are more general, the converse is not necessarily true. More- 
over, it turns out that the RG flow docs not preserve the class of models of the form 
([80|) . but takes place in the eleven-dimensional space of all vertex functions instead. 
This puts us into the favourable position of comparing the the exact flow to the trun- 
cated flow, which is given by projecting onto the space of models of the form ([50)1 after 
each stepQ In the following we present the results. 

First of all, there are three obvious fixed points, which exist both under the exact 
and the truncated flow. All three are models which can be expressed in terms of 
E- functions as in ([50)) . In these parameters, they are given by 

I) 53 BF theory: e% = e 2 Q = e 2 \ = 1 

II) Z2-BF theory: ei = 1, 620 = e 2 \ = 

III) High-temperature fixed point: e\ = e2o = e-2i = 

Point I) has only unstable directions - a slight derivation will cither lead to the model 
flowing to diverging couplings or into one of the other two fixed points. In point II) 
all couplings associated to terms containing 2g 3 are turned off. Because of (fT9")l . the 

4 Face- and edge amplitudes can in principle be absorbed into the vertex functions in this case, so 

we will fix them to equal those of the S3 BF theory, for convenience. 

5 These parameters are connected to the parameters in 1 143 II by fi\ = 1 + e\ + 2(e2o + e2i)i 02 = 
1 — ei + 2(e 20 - e 2 i), fa = 1 — a + 4(e 2 o - e 2 i), 04 = 1 + e% - (e 2 o - e 2 i). 

"There is no unique way of projecting, but different choices do not seem to yield qualitatively different 
results. 
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two remaining irrcps of S3 can be interpreted as irrcps of Z2 , so the point corresponds 
to Z2-BF theory. This point has both stable and unstable directions. In point III), 
which is completely attractive, also ls 3 is switched off. The resulting theory can be 
cither interpreted as BF theory for the trivial group, or - in the language of measure 
theory - as the Ashtekar-Lewandowski measure for S3 functions on the boundary. On 
it the kinematical and physical Hilbert space coincide, and in statistical field theory, 
it corresponds to the high temperature fixed point (HTF). 

Looking at the exact flow equations for the eleven parameters of the vertex functions, 
one can find, however, several more fixed points. These do not correspond to any 
topological theory, but can be interpreted as nontrivial. One of the points, which is 
rather close to the other three fixed points, is given, in terms of the vertex function 
couplings (i.e. the values of the {6j}-symbols) all being zero, apart from 

, S3 S3 S3 1 a 2 = { ° S3 133 ls * 

ls 3 1 s 3 ls 3 J ' 1 °s 3 1 s 3 ls 3 



e2o 



t x 1 y 










e-21 , 








tel. . 



ei 



Figure 12: Truncated RG flow of 53 model under 1—4 move. The diagram includes 
four fixed points of the flow. 



The exact recursion equations for this system are given by 

"1 > 7~~4 — ; "2 > -1 ~ 4 (81) 

1 + a\ 1 + ouy 

The obvious fixed point solution is given by a.\ = a% = 1, which corresponds to point 
II). Another nontrivial solution is e.g. given by 

IV) ai = 2.54605, a 2 = 3.31832 
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The values for a\ and a 2 can be expressed as roots of polynomials of order five, and 
we are not aware of a closed expression. The point IV) seems to have several attractive 
directions, and the correct perturbation from point I) can let the flow run into point 
IV). The point is not topological, in particular it defines a model which is not invariant 
under trivial subdivisions of bulk edges by a 2-valent vertex, unlike all the BF theories. 
Also, the resulting model is not invariant under the 2 — 3 move, as one can readily 
check. 

Since the models given by ([50)1 always lead to coefficients of the vertex functions 
satisfying a\ = e\ and a,2 = ef, one can readily check that the point IV) is not of the 
form (|50|) . It lies, however, quite close to the three-dimensional submanifold of models 
given by an _E-function. 

Under the truncated flow, in which one projects onto the space of E'-function models 
as in ([50"|) after each step, one obtains a flow in the three-dimensional space of param- 
eters (ei, e2o, e2i). One does not only recover the three obvious fixed points I), II) and 
III), but remarkably there is another fixed point which lies at roughly 

IV*) ei = 1.4021, e 20 = e 2 i = 

It can be readily checked that the vertex function derived from fixed point IV*), which 
has the only nonzero couplings a\ = e\ and a 2 = e\ , is rather close to the fixed point 
under the exact flow IV). This demonstrates that truncating the theory to models of 
the form (|80p can reveal characteristic nontrivial features of the phase space. 

6.4. Outlook 

In particular the example in the last section gives as hope that the paramctrization 
with E function presented in this work captures the relevant parameter for the rcnor- 
malization flow. More generally we have seen that there are a number of features 
similar to the behaviour of lattice gauge theories |39j . in particular the appearance 
of the strong coupling/ high temperature fixed point and the BF (low temperature) 
fixed points. (If there is a normal subgroup, a BF theory on the quotient group will 
also define a fixed point.) The BF fixed points are unstable in 2D and 3D, whereas 
the high temperature fixed points has attractive directions. It is conjectured that non- 
Abclian lattice gauge theories do not show phase transitions in 4D, hence that as in 
3D all models flow to the high temperature fixed point. 

The spin foam models provide a larger phase space and a crucial question for future 
research will be if this larger phase space allows for additional phases/ fixed points. 
This question can be investigated with spin nets models [201 115j . which are 'dimcn- 
sionally reduced spin foams'. Interestingly spin nets can be formulated with the same 
algebraic ingredients as spin foams, i.e. any spin foam model formulated with E func- 
tions leads to a corresponding spin net model. The advantage is that these models 
are already non-trivial in two dimensions which simplifies enormously coarse graining 
considerations. Additionally, it is known that 4D lattice gauge theories show a similar 
flow behaviour as the corresponding 2D (spin net / edge ) models [H] . If this similarity 
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extends to spin foam models, these investigations in 2D could already give information 
on the 4D spin foam models. 

In this work we discussed n — 1 Pachner moves, whose corresponding renormalization 
flow can be obtained exactly if vertex amplitudes are used. Classical Regge calculus is 
invariant under n—1 Pachner moves, but 4D Regge calculus is not invariant under the 
3 — 3 move [JH]- Finding theories (fixed points) invariant under the n — 1 Pachner moves 
could therefore lead to a quantum theory mimicking better the symmetry properties 
of the classical one. 

On the other hand the flow generated by the n — 1 Pachner moves is quite special. 
It leads not only to a closed, exact flow for the vertex amplitudes but also features 
invariant submanifolds, resulting from the set of vertex amplitudes that are zero if a 
given representation appears as argument. Correspondingly one can find a number of 
fixed points for this flow, which do not appear i.e. in the Migdal Kadanoff scheme 
for lattice gauge theories. The question arises, whether these fixed points have any 
significance, i.e. would flow to other fixed points, if we consider more general coarse 
grainings. 

Therefore also the other Pachner moves have to be considered, or alternatively the 
blocking of a regular lattice, which incorporates the other Pachner moves. To this 
end truncations have to be used, as otherwise the number of coupling parameters will 
grow without bound. A truncation which is suggested by this work, is given by vertex 
amplitudes that can be formulated through .E-functions. Again the main question will 
be whether this truncation captures the relevant (local) couplings. 

In this work we considered cither coarse graining procedures which only lead to local 
couplings (two dimensional theories, n — 1 Pachner moves) or a truncation of the flow to 
local couplings (Migdal Kadanoff procedure). For some questions, like the restoration 
of diffcomorphism symmetry via coarse graining [501 112] , or in order to obtain a more 
reliable phase diagram, the consideration of non-local couplings is necessary. The 
difficulty in dealing with non-local couplings is that these will proliferate and are hard 
to interpret and to control. 

One way to give an interpretation of truncations to parameter spaces including 
non-local couplings, is to introduce non-elementary building blocks which carry more 
boundary data than the building blocks one starts with [48] . The growing number 
of coupling parameters is now reflected in the growing number of boundary data, be- 
tween which the couplings can occur. These couplings can be encoded into a functional 
which acts on the boundary Hilbert space of the (non-elementary) building block in 
question. The universal boundary Hilbert space defined in section f3 . 2 1 provides in this 
sense a more general parameter space, describing also models with non-local couplings. 
(Here non-local is to be understood to include more than just the nearest neighbour 
couplings.) Again to obtain a closed coarse graining flow and not a flow which rather 
maps to larger and larger Hilbert spaces (leading to the continuum limit, i.e. a path 
integral over a region with continuum boundary data), one has to implement trun- 
cations. Such truncations can be described by maps |48| which embed the boundary 
Hilbert spaces into each other. 
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7. Summary 



In this work we proposed a new holonomy formulation for spin foam models, which is 
able to describe alj3 the current models. This class of models features, compared to 
other formulations, a minimal set of parameters, and is naturally extendable to finite 
groups and arbitrary two-complexes. We related this formulation to other holonomy 
formulations and the operator formulation of spin foams. Further applications of this 
holonomy formulation will appear in [28, 29 

The choice between different formulatiorlfl of a class of models becomes most relevant 
for coarse graining and the construction of renormalization group flows, as it defines 
the space of models in which the renormalization flow takes place (or to which it is 
projected/truncated). Here the question is still open which spin foam formulation is 
the most appropriate, at least for the four-dimensional case. 

The holonomy formulation suggest to implement coarse graining techniques from 
lattice gauge theories, and we followed this strategy for the 2D and 3D dimensional 
theories, that are in a certain sense equivalent to gauge theories. On the other hand 
spin foams have a richer structure than standard lattice gauge theories: in the E 
function formulation introduced here this can be made precise, as in spin foams E 
functions are allowed which are not class functions. The main question for future 
research will be how relevant under coarse graining the parameters are that describe 
the deviations of the E functions from a class function. This will reveal whether the 
algebraic structures specific to spin foams could lead to macroscopic phases different 
from those of standard lattice gauge theory. 



A. Non-factorizable model for a 3D dual simplicial 
lattice 

We will prove in this appendix that not all 3-valent edge C functions coming from E 
have the factorization property (|54[) . 

For a three valent edge in case of SU(2) C Spin(4) we know that every C function 
can be written in the form 

C(gi, 9 2,93) = Y, d it,jIJi,H,ji,H ^,±(51,32,33) (82) 

if 

where gi € Spin(A) and 

^ = ^^DMt\D n (g 2 )tlD n {g,)tli ]f m[m , m , (83) 

7 With the exception of the KKL model [43] , 

8 Holonomy and spin representations underly in lattice gauge theory the weak and strong coupling 
expansion respectively. Hence these different representations help to identify the appropriate 
degrees of freedom in the different regimes. For pure gravitational spin foam models such a 
distinction of regimes is much less clear, but could corresponds to a phase where geometry is 
mostly degenerate (strong coupling) and a (near) topological phase (weak coupling). 
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and Clebschcs i are normalized. 

Suppose that we are starting from the E function with some &(j+j-)h be equal 

en i\n = 1 for j ^ 1 
[3,3)0 ^ ( 84 ) 

6(1,1)1 = 1 

and otherwise 0. 

Then we can compute C -± from E functions using (|51[) as follows 



(n %tiT)k i ) bit 1 li p ut 1 > 



where fcj are the only k when e is nonzero and Pj+j-k is the projection onto SU(2) 
subrepresentation with spin k. This in fact is the product of e {j + j~)k- an< ^ a sc l uare 
of suitable contraction of Clebsches since 



< L it in p (^,.7) fci isf) = 



n 



2 



(86) 



One can check the following properties: 

1 if ji 7^ 1 and (ji, J2, J3) satisfied triangle inequalities, is an integer then 

Qiihifa,fo,h,j3 • (87) 

2 we have 

Cji ,h,h,h, i-i = > ( 88 ) 
for j'j ^ 1 because ki do not satisfy triangle inequality. 

3 We have 

( i I) . (89) 

proved by use of diagrams. 
Suppose that C-± has the factorization property then it can be written in the form 

Qj'l ,31 ,32 J2.j3.j3 = e jl.ji e h,32 e 33,33 (90) 

but from (|57|) we know that 

c w •' I) lor./ I . (91) 
By the property ([88)1 we thus obtain 

ei.i = (92) 

but then (|89[) cannot be satisfied. 

We thus proved that for this E the function C cannot be written in the factorized 
form. 
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